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Abstract 

A state vector description for relativistic resonances is derived from 
the first order pole of the j-th partial 5-matrix at the invariant square 
mass value sr = (m—iT/2) 2 in the second sheet of the Riemann energy 
surface. To associate a ket, called Gamow vector, to the pole, we use 
the generalized eigenvectors of the four-velocity operators in place of 
the customary momentum eigenkets of Wigner, and we replace the 
conventional Hilbert space assumptions for the in- and out-scattering 
states with the new hypothesis that in- and out-states are described by 
two different Hardy spaces with complementary analyticity properties. 
The Gamow vectors have the following properties: 
-They are simultaneous generalized eigenvectors of the four velocity 
operators with real eigenvalues and of the self-adjoint invariant mass 
operator M = (i-^-P^) 1 / 2 with complex eigenvalue y^s^. 
- They have a Breit- Wigner distribution in the invariant square mass 
variable s and lead to an exactly exponential law for the decay rates 
and probabilities. 
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1 Introduction 



This paper is the first of two papers which give the theoretical and mathe- 
matical foundations and the detailed derivations of certain results that were 
used in a preceding paper on the mass and width of relativistic unstable 
particles Q. In particular, it gives the construction of relativistic Gamow 
vectors which were used in |T| to describe the states of a relativistic reso- 
nance. This construction is done in complete analogy to the construction 
of the non-relativistic Gamow vectors in the Rigged Hilbert Space (RHS) 
formulation of time asymmetric quantum mechanics 0, |3|]. We shall briefly 
review some of the key properties of the non-relativistic theory in the next 
section. However, the relativistic theory does not require the non-relativistic 
theory as a backdrop and can be developed starting with the representation 
spaces of the Poincare group. 

The existence of a fundamental time asymmetry in quantum physics has 
been noted in the literature, e.g., M and ||, and the irreversible character 
of quantum mechanical decay is fairly well known ||. However an exact 
mathematical theory of decay that accommodates this property, especially in 
the relativistic domain, did not exist. Time asymmetric boundary conditions 
were incorporated in scattering theory, often unwittingly, by the heuristic 
Lippmann-Schwinger equations |7j for the in- and out- plane wave states 
\E,b k ). These were represented by Dirac kets which were supplementary 
to an otherwise time symmetric theory in Hilbert space, where the time 
evolution is necessarily reversible and given by a unitary group. In contrast, 
the time asymmetric quantum mechanics in Rigged Hilbert Spaces [0, |3| 
provides a precise definition of the Lippmann-Schwinger Dirac kets \E,b ± ) 
as functionals in a pair of RHS's of Hardy class and defines the non-relativistic 
Gamow vectors using the analytic extension of the Lippmann-Schwinger kets 
into the complex energy half-planes. 

The standard relativistic quantum theory is based on the unitary irre- 
ducible representations (UIR) of the Poincare group [| which describe 
stable elementary particles. More complicated relativistic quantum systems 
are described by direct sums (towers) or by direct products (combinations) 
of UIR's. Within this frame of unitary representations of the relativistic 
space-time symmetry group, resonances cannot be described as elementary 
autonomous physical systems. 

Our construction of the relativistic Gamow vectors as vectors associated 
to the S'-matrix pole at s = s# will lead to irreducible semigroup repre- 
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sentations of causal Poincare transformations. Like Wigner's unitary group 
representations (j, m 2 ) of space-time transformations, these semigroup rep- 
resentations are characterized by two numbers (J,Sr), the spin j and the 
invariant mass squared Sr = (Mr — zr^/2) 2 . The semigroup representations 
describe relativistic quasistable particles by relativistic Gamow vectors. Fur- 
ther, the transformation properties of these state vectors under Poincare 
transformations will show that their time evolution is exponential, and that 
the coefficient of this exponential time evolution is precisely r/2 (Tr/2 in 
the rest frame) from which we expect that the lifetime is r = Tl/Tr. Using 
Gamow vectors for the Z-boson will thus define its width as = — 2Im v /sR 



and remove the ambiguity in the definition of its mass and width [10], |TT . 

The parameters (j, Sr) that characterize the semigroup representations 
should be related to a definition of the relativistic resonance. We therefore 
start with the most widely accepted definition of a resonance as the first 
order pole of the partial S- matrix Sj(s) with angular momentum j (perhaps 
spin-parity j n ) located at s = sr (in the second sheet). We then analytically 
extend the relativistic Lippmann-Schwinger kets, \j, s, b ), from the "phys- 
ical" values So < s < oo into the complex plane and define the relativistic 
Gamow kets \ [j, Sr], b~) by the integral around the resonance pole at s = Sr. 

These relativistic Gamow kets are generalized eigenvectors of the self-adjoint 

,1/2 



mass operator (P^P 11 ) 1 with complex eigenvalue s/sr. They provide the 



state vector description of the physical entity described by the pole term in 
the relativistic S- matrix, i.e., the relativistic Breit-Wigner amplitude. 

Phenomenologically one often takes the point of view that resonances or 
decaying particles are autonomous quantum physical entities characterized 
by two real values, either the resonance mass m and resonance width T JTTJ] , 
or the mass m and lifetime r. One always assumed the relation r = h/T, 
though this could be established only in the Weisskopf-Wigner approximation 
and even the validity of the exponential law was questioned. 

Stability or the value of the lifetime does not appear to be a criterion for 
elementarity of a relativistic particle. Therefore, a mathematical framework 
capable of both uniting the notions of Breit-Wigner resonance and exponen- 
tially decaying particle and characterizing stable and quasi-stable states on 
the same footing has a strong theoretical appeal. 

Such a characterization has been accomplished for the non-relativistic 
case ||, where a decaying state has been described by a generalized eigen- 
vector of the (self adjoint, semi-bounded) Hamiltonian with a complex eigen- 
value zr = Er — iTr/2 in much the same way as stable particles are described 
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by eigenvectors of the Hamiltonian with real eigenvalues. The analogous case 
for the relativistic theory will be discussed in this paper. 



2 Key Features of Non-relativistic Gamow 
Vectors 

In the standard Hilbert space formulation of quantum mechanics, such vec- 
tors with complex eigenvalues do not exist, and for their formulation in the 
non-relativistic case one had to go to the Rigged Hilbert Space descrip- 
tion of quantum mechanics. This was not a revolutionary step since the 
mathematical definition of Dirac kets already required the Rigged Hilbert 
Space. The generalized eigenvectors with complex eigenvalues, ip G = \Er — 
iTji/2~) V271T, which we called Gamow vectors because of their property 
(|2.3| ) below, were then obtain from the pole term of the S'-matrix. The real 
parameters Er and Tr are respectively interpreted as resonance energy and 
resonance width, for reasons that will become clear from their property ( ]2.3| ) 
below. 

Like Dirac kets, Gamow vectors are functionals of a Rigged Hilbert Space 
(Appendix A): 

$+ C H C : ip G G (2.1) 

and the mathematical meaning of the eigenvalue equation H x \z R ) = zr \z r ) 
is: 

{H^\z r ) = {,Ij\H X \zr) = zr^\z~ r ) for all ^ G $+ . (2.2) 

The conjugate operator H x of the Hamiltonian H is defined, by the first 
equality in ( [2.2|) , as the extension of the Hilbert space adjoint operator 
of H to the space of functionals Q 

The non-relativistic Gamow vectors have the following properties: 

1. They have an asymmetric time evolution and obey an exponential law: 

ifF(t) = e+ iHXt \E R - ir/2") = e- lE ^e- vt ' 2 \E R - iT/2~), 

only fort>0. (2.3) 



1 For (essentially) self-adjoint H , H> is equal to (the closure of) H; but we shall use the 



definition (2.2) also for unitary operator U where U x is the extension of and not U. 
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There is another Gamow vector ip G = \E R + iY /2 + ) e $*, and a 
semigroup eZ tHXt for t < in another Rigged Hilbert Space $_ C 7i C 
(with the same 7Y) with the asymmetric evolution 

= e -iH*t\ ER + ir / 2 +) = e - iE ^e rt / 2 \E R + iF/2+), 

only fort<0. (2.4) 



2. The ^ G (V ,G ) is derived as a functional associated to the resonance pole 
term at z R = E R — iY/2 (at z* R = E R + iY/2) in the second sheet of the 
analytically continued S'-matrix. 

3. The Gamow vectors have a Breit-Wigner energy distribution 




<~ g| ^ ) = 'Vfr g-(g fl -g/2r -°°" <g< °°- (2 ' 5) 

where — oo/j means that it extends to — oo on the second sheet of the S- 
matrix (whereas the standard Breit-Wigner terminates at the threshold 
E = 0). 

4. The decay probability P(t) = TriA v \ip G (t))(ip G (t)\) of ip G (t), t > 0, 
into the final non-interacting decay products of the channel described 
by can be calculated as a function of time. From this the decay rate 
R(t) = is obtained by differentiation as R(t) = r^)e~ rt where 
is the partial width for the decay channel rj. And it leads to an exact 
Golden rule for Y^ which in the Born approximation (ip G — > f D , an 
eigenvector of H Q = H — V; Y/E R — * 0; E R — > E ) goes into Fermi's 
Golden rule No. 2 of Dirac. As a consequence of the exponential law for 
R(t), the lifetime of the state described by ip G (t) is given precisely by 
r = h/Y. 

We want to generalize these non-relativistic Gamow vectors to the rela- 
tivistic case and construct representations of the Poincare group V that de- 
scribe relativistic resonances and decaying states very much in the same way 
as unitary irreducible representations (UIR) of V describe stable relativistic 
particles. In order to obtain a state vector description for unstable particles, 
one has to start from the space of decay products. For a relativistic unsta- 
ble particle decaying into a two-particle system with masses mi, m 2 , and 
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spins si, s 2 , the (asymptotically free) decay product space is the direct prod- 
uct space of irreducible representation spaces of the two particles 7i(mj,Sj), 
i — 1, 2: 7ii2 = 7i(mi,si) ® 7i(m 2 , S2). A decaying state vector can then 
be obtained, in analogy to the non-relativistic case, by extending the basis 
vectors of TC12 to complex energy, in a manner which is consistent with the 
dynamical description of unstable particles as poles of the analytically con- 
tinued S'-matrix. Since resonances occur in one-particular partial-wave j^, 
(which designates the spin (total angular momentum in the rest frame) and 
parity of the unstable particle), we need to use the basis vectors of Hu which 
are diagonal in total angular momentum in order to obtain the relativistic 
Gamow vectors. These angular momentum basis vectors are obtained by 
the reduction of the direct product TC(mi, si) ® TL{m 2l s 2 ) into a continuous 
direct sum of irreducible representation spaces 

H 12 = J2 dfM(5)Hl(s,j), (2.6) 

• J(m, 1 +m 2 ) 2 

where s is the invariant mass square (Mandelstam parameter) for the two- 
particle system s = (pi + P2) 2 , and r\ and n are degeneracy and particle 
species labels respectively |13| . As discussed in Sections [T2]and|4] below, 
the partial S'-matrix is a function of s which is analytic except for poles and 
branch cuts. The complex poles at s = s# are the ones associated with the 
unstable particles, and to these values the analytic extension in the s- variable 
will be performed. 

The extension of the invariant mass square s to a complex value leads 
necessarily to complex momenta, since s = E 2 — p 2 . As Lorentz transfor- 
mations intermingle energy and momenta, this in general leads to complex 
momentum representations of the Poincare group. To obtain a description 
of unstable particles with minimal modifications to the stable particle case, 
and to restrict the set of complicated complex momentum representations of 
V jl^l ) we will consider representations of V where complexness is due only to 
the complexness of the mass square s. These representations are "minimally 
complex", in the sense that, while the invariant mass squared s is complex, 
the 4-velocities = p M / ^/s remain real. This construction was motivated 



by a remark by D. Zwanziger [Tj| and is based on the fact that the velocity 
kets provide as valid a basis for the UIR's of the Poincare group as Wigner's 
momentum kets. Moreover, the 4-velocity eigenvectors are often more useful 
for physical reasoning, because 4-velocities seem to fulfill to a rather good 
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approximation "velocity super-selection rules" which the momenta do not 
The reduction of TI12 in the velocity-angular momentum basis \pj3[sj]r], n) 



0. 



was performed in |13[ as a preparation for this paper. The discussion there 
extends to an unstable particle decaying into more than two particles, but 
the reduction of the decay product space into a direct sum of UIR spaces 
becomes much more tedious than for the two-particle case. Hence, in gen- 
eral, with the use of velocity-angular momentum kets, the relativistic Gamow 
vectors will be defined not as momentum eigenvectors but as 4- velocity eigen- 
vectors in the direct product space of UIR's for the decay products of the 
resonance R. This definition requires the extension of the invariant square 
mass s = (pi + p 2 + ■ ■ ■ ) 2 , where pi, p 2 • • • are the momenta of the decay 
products of R, into the complex values. In analogy to the non-relativistic 
case, we define the resonances R by poles of the analytically continued rel- 
ativistic partial 5*-matrix with total angular momentum j = jr and pole 
position s = Sr = (Mr — iT R /2) 2 . As discussed in Section |3j] below, there 



appear to be phenomenological and theoretical arguments in favor of this 
definition. This means that the relativistic particle is characterized by the 
value of its spin jr and its mass wr = v /sr = Mr — iT R /2, where T R is 
zero for stable particles, but wr is complex if the particle is unstable. Just 
as the mass and spin of a stable particle permit its association to an irre- 
ducible unitary representation of the Poincare group characterized by these 
values, it will be shown in a forthcoming paper that the complex mass Wr 
and spin jr concatenate an unstable particle to an irreducible representa- 
tion of the Poincare group, only these irreducible representations are not 
unitary. In fact, it turns out that they are irreducible representations of the 
causal Poincare semigroup, defined as the semi-direct product of the group 
of homogeneous Lorentz transformations with the semigroup of space-time 
translations into the forward light cone. The fundamental mathematical ob- 
ject needed in constructing these semigroup representations is the relativistic 
Gamow vector. The main technical result of this paper is its derivation, while 
the discussion of the semigroup representations is deferred to a later paper. 

3 Relativistic Gamow Vectors From Poles of 
the Relativistic S'-matrix 
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3.1 Introduction 



As discussed above, our aim is to obtain the relativistic Gamow vectors 
from the pole term of the relativistic S'-matrix in complete analogy to the 
way the non-relativistic Gamow vectors were obtained |T7|]. For the sake of 
definiteness, we discuss here the case of resonance formation in an elastic 
scattering process ah — > R —>■ ah for which one may consider 7r + 7r~ — ► p — > 
18|| or e + e~ — > Z — > e + e~ [jnj as examples. 



7T 1 7T 



In the absence of a vector space description of a resonance, we shall 
also in the relativistic theory define the unstable particle by the pole of 
the analytically continued partial S-matrix with angular momentum j = jr 
at the value s = Sr ; = (Mr — f Tr) 2 of the invariant mass square variable 
(Mandelstam variable) s |1{J. This means that the mass Mr, width Tr (which 
will be shown to be connected to the lifetime by tr = H/Tr), and spin jr 
are the intrinsic properties that define a quasistable relativistic particle Q. As 
will be shown below in Section |5|, the j-th partial S'-matrix can be separated 
into a background term and a pole term which is given by the relativistic 
Breit-Wigner amplitude 

a,- = / n R ^ . — with — oo/j < s < oo , (3.1) 

m s-{M R -\Y R y ii- - , 

where the negative values of s are on the second sheet. 



2 The more conventional definition, used in fllTf , of resonance mass and width is M p , 
T p , which are defined in terms of Mr, Tr by 

/ i / p \ ^ 



s R ee Ml - iM p T p = M 2 R ( 1 - - ^ ) ) - iM R T R (|yb) 



or M p = Mr^J 1 - | (i^) 2 and T p - r R ( 1 - \ (^)'J • M p 2 is the peak position 

of the relativistic Breit-Wigner probability |a_, 3 (s)| 2 . The exponential time evolution of the 
Gamow vectors which will be derived in the forthcoming sequel to this paper, shows that 
r r, a nd not T p , is the inverse lifetime. Since only Tr (not Tz or T p ) of the parameteriza- 
tion ( j3.l| ) precisely fulfills the identity Tr — K/t, we call Tr the "width" of the relativistic 

resonance. For the p-meson (m^) ~ 0-03 and for the Z-boson (j^j ~ 7 x 10 -4 . 

Usually the difference between M p and Mr is an order of magnitude smaller and/or well 
within the experimental error for the experimental mass; however for the p-meson and in 
particular for the Z-boson this difference is significant, see p0|. 



-1/2 
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However we want to mention that the S-matrix pole definition with the 
parameterization by Mr, Tr is neither the only nor the universally accepted 
definition of the parameters mass m and width T of a relativistic resonance. 
The particle data table uses M p and T p as mass and width 2 . And for 
precision data as in the line shape analysis for the Z-boson [fTofl , one prefers 
an energy dependent width replacing M p T p — > v /ir(s), cf. equation (35.53) 



of [11]. This gives still another parameterization which in the case of the 



Z-boson is called Mz, T 'z ||10|| . In the line-shape analysis for the Z-boson, 
these three parameterizations lead to three different values for the Z-boson 
mass and width GJJ which fit the experimental data equally well. Therefore 



the line-shape data do not provide a phenomenological clue as to how the 
fundamental parameters m, T of the Z-boson should be defined EH]. In 



the case of the p-meson, the definition of the resonance parameters given by 
( |3.1j ) leads to better, process- independent fits of the parameters KJJ . Thus, 



there seems to be a slight phenomenological preference for the S'-matrix pole 



definition ( |3.1|) of a relativistic resonance. The parameterizations in terms of 
Mr, Tr or M p T p of the S'-matrix pole position sr are phenomenologically 
equivalent. One of the results of our theory will be that only Tr and not 
Tp (or Tz) will be connected to the lifetime by Tr = H/tr, which is the 
theoretical argument in favor of the parameterization sr = (Mr — iTR/2) 2 . 

In order to make the analytic continuation in the partial S-matrix with 
angular momentum j, we need the angular momentum basis vectors for the 
scattering states. Thus we replace the non-relativistic angular momentum 
basis vectors in the derivation for the non-relativistic Gamow vectors with: 

\E) = \Elh V }^\pj 3 [w = V~s,j}} (3.2) 
| E ± ) ee | Ellrf > -> | pj 3 [w = V~s,j} ± ). &) 

The two sets of bases are related by the Moeller wave operators 

I PJs[w = y/s, j] ± ) = Q ± \pj 3 [w = y/s,j] ) . 



We shall analytically extend the Dirac kets (|3.2p ) in the variable s from 
the physical values to the complex values on the Riemann surface of the S'- 
matrix. The analytic extension of \p,j 3 , [s, at the pole position Sr yields 
the relativistic Gamow vector. 



9 



3.2 The S-matrix 3 



In a scattering experiment, the experimentalist prepares an initial state <f) m 0, 
describing the non-interacting projectile and target beams, at t — > — oo. 
Then, later, at if — ► +oo the experimentalist measures or registers an ob- 
servable 4 \%l)° ut )(ip out \. It is assumed that the time translation generator H 
can be divided into two terms, the "free-particle" Hamiltonian K(= P® + P 2 ° 
at rest) and an interaction part V: 

H = K + V, 

where the split of H into K and V will be different if different in- and out- 
particles are involved. 

The state vectors 4> m (t) = e~ lKt (j) m and the observable vectors ip out (t') = 
e~ ip out evolve in time according to the free Hamiltonian K. When the 
beams reach the interaction region, the free in-state vector m turns into the 
exact state vector <p + whose time evolution is given by the exact Hamiltonian 
H = K + V: 

tt + <p m (t) = 4> + (t) = e - iHt 4> + = Q-4> out (t) . (3.3) 

Here t is the proper time in the center-of-mass of the projectile and target. 
This vector <p + leaves the interaction region and becomes the well determined 
state <f) out . The state vector <p out is determined from <p m by the dynamics of 
the scattering process: 

^out = s< pn^ g = ( 3 _ 4 ) 

The state </> m and thus <ft + and also <p out are determined by the preparation 
apparatus (the accelerator). 

A scattering experiment consists of a preparation apparatus and a regis- 
tration apparatus (detector). 

3 Here we closely follow chapter 3 of |j in order to both display the analogy (and compa- 
rability) and expose the differences between our development and the standard views in rel- 
ativistic quantum theory. Our notation transcribes into that of || as {<f> m / out , ip out } — > $ g 
and {<p + ', Tp~} — > . In Q, the multi-particle basis vectors are also denoted by where 
a = {pi<Jini 1 p2<J2n2 1 ■ ••}, cr is the third component of the spin, and n is the species label. 

4 In realistic experiments the states are not pure but mixtures W m = J2 w a\4'a l )(4 >t a\ 
and the observables are not given by projection operators but by A ouf = 

EW)\1>b*)W*\ 
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The registration apparatus registers an observable \ip out )(if) out \ outside the 
interaction region. This observable vector ip out comes from a vector 

rj,- = VL~^ out (3.5) 

in the interaction region. The observable ip out is of course not the same as 
the state <f) out , since (p out , and thus + and <p m , is defined by the accelerator 
and ip out , and thus ip~ , is defined by the detector. Hence he set of vectors 

may be, and in our case is, distinct from the set of vectors (f) + . 

The probability amplitude to register the observable \ip out ) (ip out \ in the 



state <j) out , which-by simple calculations using (3^3) and ( |3.4|) -is the same as 
the probability amplitude for the observable in the state + , is the 

scalar product 

(ij out ,0 out ) = (i, ou \s<p m ) = (n-ij out ,n + m ) = (i>-,<f> + ) . 

The fundamental probabilities of quantum theory (Born probability) are the 
probabilities for the observables \ip~)(ip~\ in the state (or for the 

observable vector ip~ in the state vector + ) given by 

Tr(|^-)(^|0 + )(0 + |) = |(^,0 + )| 2 . 
The vectors 0+ are called the in-states and the vectors ipg are called 
the out-states and the array of complex amplitudes (ipg, 4>a) ls called the S- 
matrix. The labels a and f3 stand for a whole collection of discrete quantum 
numbers. The S'-matrix is also defined when a and/or (3 are continuous labels 
(of basis vectors for the space of the {if)~} and the space of the {</> + }) only 
then the S'-matrix does not represent probability but a probability density. 
"It should be stressed that in-states + and out-states do not inhabit 
two different Hilbert spaces" || but they do inhabit two different (dense) 
subspaces of the same Hilbert space H; these two dense subspaces we call 
$_ = {0 + } and $ + = {f>~}- This is the only new hypothesis by which 
our time asymmetric quantum theory in Rigged Hilbert Space differs from 
the standard Hilbert space quantum theory. Thus the Rigged Hilbert Space 
theory distinguishes meticulously between prepared states (in-states) {4> + } 
and observables (out-states) by means of two different Rigged Hilbert 

Spaces of Hardy class: 

+ G C Tt C for the prepared in-states defined by the 

preparation apparatus (accelerator) _ ^ 

ip~ G $+ C Ti. C $^ for the registered out-states (observables) 

defined by the detector 



11 



In contrast, the conventional scattering theory assumes $_ = $+(= TC)- As 
a consequence of this new hypothesis, an arbitrary in-state <p + G $_ can- 
not be expanded as a sum of out-states ip~ G However (for instance 
in the non-relativistic theory), any in-state <fr + can be expanded with re- 
spect to a generalized basis system \E~) G of the out-states ip~ (ip~ = 
J dE\E~)(~ E | using the 5-matrix : 

/"OO />oo 

</>+=/ dE\E-)(~E\<p + ) = dE\E~)S(E + iO)( + E \ </> + ) . 
Jo Jo 

In the non-relativistic theory the boundary conditions required of the 
spaces $_ and $+ lead to distinct analyticity properties for the energy wave 
functions ( + E | <p + ) of <p + G $_ and the energy wave functions (~E \ of 
•ty- G $+ H- To wit 

G $+ if and only if = e5fl H 2 + \ R+ (3.7a) 

0+ G $_ if and only if = {E\(f) in } eSf]H 2 _ | R+ (3.7b) 

where: 5 denotes the Schwartz space; H±, the Hardy class functions (Ap- 
pendix 0) which are boundary values of functions that are analytic in the 

open na lf plane of the complex energy plane; and |r + , the restriction 

to the physical values of the energy, IR + = [0, oo). As the lower half complex 
energy plane we take the second Riemann sheet of the S"-matrix reached from 
the physical values ( + E + iO \ <f) + ), E G K+ by burrowing down through the 
cut along the positive real axis R + . As the upper half complex energy plane 
we take the second Riemann sheet of the S'-matrix reached from the values 
(~E — iO | E G M + on the physical sheet by burrowing up through the 
cut along R + . In the non-relativistic case ( |3.7b| ) and ( |3.7a| ) could be obtained 



from a mathematical formulation of causality (using the Paley- Wiener the- 
orem) p2"| . An important consequence of the distinction between $ + and 
$_ obtained in ( |3.7| ) is that the time evolution operator U(t) = e lHt , which 



is represented by a group on the Hilbert space Tt = L (M + ), splits into two 



semigroups when restricted to <&± [17]. This means that U(t)\<s, + is a contin- 



uous operator on $ + only for t > and Z7(t)|$_ is a continuous operator on 
$_ only for t < 0. It is this semigroup property that is at the origin of the 
time asymmetric evolution equations (|2.3|) and ( |2.4j) . 

In the relativistic theory, in particular for the elastic scattering of two 
particles, the Hilbert space TC of (|3.6|) is the direct integral (|2.6|) of irreducible 
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representation spaces 7"^(s, j) labeled by the Mandelstam variable s = (pi + 
P2) 2 and the total angular momentum j of the two particles. For the measure 
dfi(s) in ( |2.6| ) we take the Lebesgue measure ds. In what follows, we will 
restrict our discussion to the elastic scattering of two spinless particles with 
the same mass m. We use as a basis for $-t the angular momentum velocity 
basis vectors, {|pj3[sj] T )} [[CJ. With these basis vectors, the Dirac basis 
vector expansion for an in-state (of two equal mass spinless particles) <p + G 
$_ is: 

^ = E r ds l SlPJ3[sj] + )(pj 3 [sj] + |0 + ) , (3.8) 
333 j4m2 J P 

and for an out-observable ip~ G $ + is: 

poo f* j3 a 

^~ = E / rfs / ^||Pi3[sj]-}(pj 3 [sjriV-) • (3.9) 

„•• Jim 2 J 



In (|3.8| ) and (|3T9|), (pj3[sj] + |0 + ) and (pj3[sj] j-0 ) are the wave functions 



of and ip~ along the velocity basis vectors. We shall make the hypoth- 
esis that these wave functions have the same analyticity properties in the 
invariant mass squared s = (E cm ) 2 as the energy wave functions in ( |3.7|) . 
However, in the relativistic case, due to the mathematical requirement of 
the invariance of the subspaces Q T under the action of the generators of the 
Poincare group, a closed subspace S of the Schwartz space S of (|3.7|) has to 



be considered. The subspace S, constructed in [E3J, is the space of Schwartz 
functions which vanish at zero faster than any polynomial. This requirement 
also assures that the zero mass states do not contribute to the Gamow vec- 
tor (see (|5.17|) ). This avoids the difficulty that the four velocity operators, 
which is centrally significant to our construction of Gamow vectors, cannot 
be meaningfully defined in the zero-mass case in any obvious way. The fea- 
tures of the space S which are needed for the construction of the relativistic 
Gamow vectors are as follows 1 12311: 



Proposition 3.1. The triplets 

snn 2 T W H) c l 2 (m so ) c (<Snw|| Rso ) x (3.10) 

form a pair of Rigged Hilbert Spaces. 
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In ( B.lUj ), M so is the set of physical values of the Mandelstam variable s for 
the scattering process, M so = [(mi + m 2 ) 2 , oo). 

Proposition 3.2. The space S is endowed with a nuclear Frechet topology 
such that multiplication by st ? 

s% ■. snnl^snnl, n = 1,2,3,- •• 

is a continuous linear operator in the topology of S . 

Thus the relativistic characterization of $± analogous to ( |3.7a|) and (|3.7b|) 



is: 



e $ + if and only if (pj 3 [sj]~ \^~) G5fl x 5 ( R ) ( 3 - lla ) 

0+G$_ ifandonlyif (pj 3 [sj] + |0 + ) G5n^ 2 | Rso x5(M 3 ); (3.11b) 

where R 3 is the space of components of the 4- velocity and the Hilbert space 
7i of (|3.6|) is realized by the function space 



rfs)xL 2 (M 3 ,0) • (3.12) 



In ( |3.11| ) as in ( |3. 7] ), 7i+ means the functions of Hardy class analytic in the 
upper half of the second sheet of the s-plane and 7i 2 _ means the functions 
of Hardy class analytic in its lower half. Specifically, the physical values 
( + s — iO | cf) + ) are the boundary values of functions analytic in the lower half 
of the second sheet and the (~s+iO | ip~) are the boundary values of functions 
analytic in the upper half of the second sheet. These analyticity properties 
on the second sheet of the complex s-Riemann surface will turn out to be 
important because the (pairs of) resonance poles of the S'-matrix are located 
on the second Riemann sheet. 

By virtue of Proposition |^, the total momentum operators P^ = P\^ + 



P2fi and the invariant mass square operator M 2 = P^P M are r$ ± -continuous 
operators; hence their conjugates P* and M 2 *, are well defined on 
This can be seen by considering the realization, for instance, of the vectors 
P^- and M 2 ip~: 

(P^-\pj 3 [sj]-) = ^-\P^\pj 3 [sj]-) = V^Mi>-\ph[sj}-) , (3.13a) 
(M 2 i>-\pj 3 [sj}-) = ^-\M 2X \pj 3 [sj]-) = s{i/>-\pj 3 [sj]-) . (3.13b) 



5 defined by the first equality in (3.13) 
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According to Proposition |3.2| and the definition of the wave functions Ipjafsj] - ) 
given in (|3.11a|) , the multiplication operators by \/sp^ and s which appear in 
the right hand side of ( |3.13a| ) and ( |3.13b| ) are r,j, + -continuous. Consequently, 
P^ and M 2 are r$ + -continuous operators, and the conjugate operators M 2 * 
and P* that appear in ( |3.13| ) are everywhere defined, weak* -continuous op- 



erators on $*. Hence, (|3.13a|) and ( |3.13b|) define the functionals |pj3[sj]~) 
as generalized eigenvectors of P M and M 2 . The same discussion applies for 
the space Summarizing 

P^ : $-t — > $± is r$ ± -continuous (3-14) 

and 

P^\P3 3 [sjD = VwnlPMsj]*) , (3.15a) 



M 2X \pj 3 [sjF) = s|pj 3 [sj] =F > 



(3.15b) 



We can re-express the generalized eigenvalues of the momentum operator in 



terms of the three velocity v by noting that p = "yv 



\/l—v 



2- 



and 1 + p 



l-v 2 



7 2 . Hence, the eigenvalues in (|3.15a|) can be rewritten as 



H x \pj 3 [sj} T ) = jV^\m[sjF) , 

P X \P33[sj} T } = lVsv\p] 3 [s]] T ) . 



For the branch of y/s in (|3.13|) , (|3.15|) and ( |3.16|) , we choose 
—7i < Arg s < 7r . 



(3.16) 



(3.17) 



This choice of branch, even though irrelevant for the physical values of s, 
will be needed since we will analytically continue the kets |pj3[sj] T ) to the 
unphysical second Riemann sheet as described in Section [|. 
We shall now consider the 5*-matrix element 

(^s<n = or,0 + ) 

333 



3% 



< PJs [sj] I S | p% [s'f] ) (+p% [s'f] I <P + > (3-18) 
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where we insert into (ip , + ) a complete system of basis vectors [] and use 

(-pjsM I P'Ms'jT) = ( n-| pj 3 [sj] ),n + \p%[s'f} > ) 

= (ph[sj)\n-^ + \p'j' 3 [s'j')) (3.19) 

= (Pj3[sj]\S\p'f 3 [s'f]) 

Using the invariance of the S operator with respect to space time translations 
[P„S] = (3.20) 
we obtain 

(PjsMvlsip'j'As'j'W) 

= S(p - p')6(pa - P^dmMv I s I P'Ms'fW » (3.21) 
where (( \S\ )) is a reduced S-matrix element. This we can also write as 

( PJsMv I S I P'Ms'j'W ) = 2E(p)5(p - p')5(s - s') 

((ph[sj]v\S\p%[s'j'W)) (3.22) 

where (( \S\ )) is another reduced matrix element defined by ( |3.22|) . In 
( |3.21| ) and ( p.22[ ) we include explicitly the degeneracy quantum number rj for 



purposes of clarity and completion, but we will omit it below for the sake of 
notational convenience. The form ( |3.22| ) follows from ( |3.2lD by the defining 



identities p = ^j=, p° = From the invariance of the S'-operator with re- 
spect to Lorentz transformations, in particular from W (L~ 1 (p))SU (L _1 (p)) = 
S [] it follows that the reduced matrix element is independent of p (i.e., it is 
the same for all p as in the center of mass frame p = 0). Invariance with 
respect to rotations in the center of mass frame shows then by analogous 
arguments for the discrete quantum numbers j'3 and j that the reduced ma- 
trix element is proportional to 8j 3 j' 3 8jj' and independent of j 3 . Since Poincare 
transformations do not change the Poincare invariants s and j, the reduced 
matrix element can still depend upon s and j. Thus we have 

( Pjs[sj}r] I S I p'Ms'j'W ) = 2E(p)5(p - p')5(s - 

(i7||S,(s)||i/> (3.23) 



6 We ignore the possible existence of bound states of H of which there are usually none; 
certainly not for the ir + ir~ system of tt + it~ — > p° — > tt + it~ . 
L(p) is the rotation free boost. 
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If there are no degeneracy quantum numbers, i.e., if we ignore the particle 
species label and channel numbers and restrict ourselves to the case without 
spins (like for the tt + tc~ system), then the reduced matrix element can be 
written as 

( V \\S j (s)\\rf)=S j (s) (3.24) 

where j is the orbital angular momentum in the center of mass. We insert 
( |3.23| ) and fl3.24| ) into ( |3.18| ) to obtain for the S'-matrix element 



(*-.* + ) = E f *£ / ^<*-|£Wfl)|0»O 

S j (s)( + Oj a js\rf{L(p))\<l> + ) (3.25) 

Since resonances come in one partial wave with definite value of angular 
momentum, we consider only the one term in the sum over j with j = Jr{— 
1~ for 7r + 7r~ — > p° — > 7r + 7r~), i.e., we restrict ourselves to the subspace with 
j = Jr {s = 0, I = j, n = n p , n n7T ). This means that we consider only 
the term with j = jn in the sum on the right hand side of ( p. 251 ) and call 



Sj R (s) = S(s). To simplify the equations, we also consider ( |3.25|) first for 
fixed values of j 3 and p, say p = and j 3 = 0. We define 

| j 3 , s-> = (ip- | U{L(p)) | j 3 , j = ]rs~) 

= (s,h\^-{p)) esnn 2 _ k (3.26) 

(+0 j 3 , j R s | U\L(p)) = ( + j 3 ,s | + (p)> G 5 n ft 2 | Kso 
and write (I3.25J) 

^|fe(p),ft(p)) (3-27) 



J3 

where 



(p), 0+(p)) = J dsty- | E7(L(p)) | 0j 3 , jflS-)S'(s) 

( + Oj3,^s|^( L (p))|0+) (3.28) 

= J ds(i;-(p)\j 3 s-)S(s)( + j 3 s\ ( p + (p)) . 
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For the fixed values p = 0, j 3 = we suppress the labels p and j 3 : 

(^-,0+) = {^- =0 (p = O),0+ =o (p = 0)) (3.29) 



and write 



V>~,tf» + )=/' c/s(^"|s-)S(s)(+s|0+). (3.30) 



(rai+m2) 2 

The quantities ( , _ , + ) are not matrix elements but matrix element densities 
in the sense that they have to be calculated for all 4-velocities p and all j% 
and then integrated and summed using ( p. 27] ). Considered as functions of 



the invariant mass-squared s, we make the assumption fl3.11|) , i.e., we assume 
that the wave functions fulfill 

(+s|0+)e 5n^ 2 k o (331) 
(s | if;") = | s-) e 5 n n 2 + \ Rso , i.e., \s-)esnn 2 _ | Kso . 

Then ( |3.30| ) is very much the same as the corresponding expression for the 



or 



S-matrix element in the non-relativistic case (e.g., equation (3.9) in [22 
equation (5.4) in JJ) except that in place of the non-relativistic energy E 
with lower bound we have here the center of mass energy squared s = 
E^ m , with lower bound (mi + r/22) 2 . The reason for this choice is that s 
is the variable of the relativistic partial S'-matrix Sj(s) which is continued 
to complex values, and the Riemann energy surface of Sj(s) has features 
similar to the Riemann surface of non-relativistic partial S'-matrix Si(E). In 
particular, in the relativistic S-matrix theory Sj(s) is analytic on the first 
"physical" sheet except for cuts along the real line (and bound state poles 
if any exist), and resonances are defined by (pairs of) poles at the complex 
values s = = (Mr ± ^Tr) in the second sheet or another "unphysical" 
sheet. 



4 Analytic Properties of £(s) 

We start with the hypothesis that unstable particles and resonances are asso- 
ciated with poles of the relativistic partial S-matrix (rj || Sj(s) || rj') = S(s) on 
the "unphysical" Riemann sheet. The partial S-matrix S(s) for physical val- 
ues of the center of mass energy s > (mi + m.2) 2 is the boundary value to the 
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real axis s + ie , e — > + of a function in the complex s-plane which is analytic 
except for cuts along the real axis and possibly bound state (stable particle) 
poles Pi at s = m 2 P . < (mi + mi) 2 ■ The right hand cut (a consequence of 
unitarity) starts at the elastic scattering threshold s m i n = (mi + 1712) 2 . This 
cut is two-sheeted JL8|, i.e., of the square-root type. The elastic scattering 



matrix element, rj' = 77, has further branch points at each energy corre- 
sponding to a threshold for a newly allowed physical process (e.g., possibly for 
7r + 7r~ — > p — > 7r + 7r~7r + 7r~ at Sth = (2mi + 2m2) 2 ). To make S(s) single-valued 
on a Riemann surface, cuts start at these branch-points and are drawn along 
the real axis, cf. Figure [I]. These branch-points are called normal thresholds 
and start at energies at which production of other particles (inelastic pro- 
cesses) is possible. If one does not cross a cut, one stays on the "physical 
sheet". The first normal threshold s = (mi + m<i) 2 is the least s at which a 
two particle state can exist. To reach other sheets of the Riemann surface of 
S(s) one burrows through one or several branch cuts. These sheets are called 
the "unphysical" sheets. The second sheet is reached from the physical value 
s + ie , s > (mi + m.2) 2 by burrowing down through the normal threshold 
between s^i = (mi + m 2 ) 2 and s^ 2 . The resonance poles P of the rela- 
tivistic elastic scattering S-matrix element is located on the second sheet at 
Sr = (Mr — f r fl ) 2 . This hypothesis is suggested by the idea that an unstable 
particle may be connected to a stable particle by letting the parameters that 
measure the strength of the force (between the decay products) vary contin- 
uously until that force becomes so strong that the unstable particle becomes 
a stable particle, since a stable particle corresponds to a pole on the physical 
sheet at the real value s = mf>. below the threshold (mi +m2) 2 ; in the transi- 
tion to instability the pole must pass round the (mi + 1712) 2 branch point and 
through the cut. Hermitian analyticity (symmetry relation of the S'-matrix 
S(s — ie) = S*(s + ie)) implies that when the pole position is complex there 
must be a pole P' at the complex conjugate position = (M R + jT R ) 2 on 
the unphysical sheet reached by burrowing through the cut from the lower 
half plane of the physical sheet of Figure |l|. Thus a scattering resonance is 
defined by a pair of poles on the second sheet of the analytically continued 
S'-matrix located at positions that are complex conjugates of each other. The 
pole P 1 corresponds to the time-reversed situation which we do not want to 
discuss here. (There may exist other resonance poles located at the same or 
other physical sheets, but we will mainly be concerned here with one pair 
of poles in the second sheet). There may also be higher order poles |24] or 



branch-points on the unphysical sheet but we do not venture here into these 
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complications; for the sake of simplicity, we shall restrict ourselves to the 
elastic process tt + tt~ — > 7r + 7r~. and assume further that there are no other 
resonance poles in the same channel (ignoring a possible p(1450)). 



5 Analytic Extension to the Resonance Pole 

We want to analytically continue the matrix element density (ipj^p), 0l(?>)) 
(|3.28| ) so as to encompass the resonance pole contribution to the scattering 
amplitude. As a specific example, we consider resonance formation in elastic 
scattering, e.g., 



p° -> TT+TT . (5.1) 



We will split (see (|5.8|) below) (ipj 3 (p), <j>j 3 {p)) into a resonance term and a 
background term. The resonance term yields the state vector description 
of the unstable particle, which in this case is the p°-meson. The procedure 
depicted below can be generalized to processes other than (|5.1|) , with the 
background term being determined by the system of branch cuts specific to 
the particular process considered. 

For ( |5.1| ), + is determined from m Q3.3| ), the incoming two pion state 
prepared by the preparation apparatus, and ip~ is determined from ip out 
the registered two pion state. Since the p meson has spin 1, in (|3.28|) 



we have ]r — j p = 1- The path of integration in fl3.28|) extends along the 



positive real axis, just above the normal threshold cuts mentioned in Section |] 
from 4m 2 to oo, where m is the mass of 7r + . As explained in Section |], 
if we assume that the least energetic channel following 7t + tt~ — > tt + tt~ is 
7r + 7r~ — > 7r + 7r~7r + 7r~, then ^(s) has a square root branch cut between 4m 2 
and 16m 2 . Thus the physical sheet is connected to the second Riemann sheet 
between 4m 2 and 16m 2 and S'(s) can be analytically continued everywhere 
in the second Riemann sheet jn|, except for the resonance pole of the p° 
meson. So, to analytically extend ( |3.28| ) we deform the contour of integration 
between 4m 2 and 16m 2 in the first sheet into an infinite semi-circle in the 
lower half plane of the second sheet (cf. Figure ^|), taking into account the 
p-pole of S(s). This is the only singularity of the integrand in ( |3.28| ) since 



the wave functions (ip (p) \ j 3 s ) and ( + jss \ (f> + (p)), according to (|3.11 



are 



Hardy class from below, hence analytic for Ims < 0. Explicitly, the contour 
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deformation of ( |3.28| ) yields, according to Figure 



{i>J 3 {P)At{p)) 



33 1 



ds(^-(p)|j3S-)^(s)( + J 3 s|0 + (p)) 



4m 2 



4m 2 



ds(ij-(p)\j 3 s-)S H (s - ie)( + j 3 s\4> + (p)) 



+ I ds(r(P)\j3S-)S II (s)( + j 3 s\ ( f> + (p)) 

J Cno 



+ yds^-(p)|j3S-)5 // (s)( + j3S|0 + (p)) 



+ 



16m 



ds^-(p)|j 3 s->5"(s - ^e)( + J3s|0 + (p)) 



+ 



16m 2 



^-(p)|j3S-)5(s + ^)( + j 3 s|0 + (p)) . (5.2) 



In the above expression, the terms between the brackets [ ] result from the 

analytic continuation of f 4 ™ ds ■ ■ ■ into the second Riemann sheet through 
the square root cut between 4m 2 and 16m 2 , S 11 is the S-matrix in the sec- 
ond sheet, Cqo refers to the infinite semi-circle in the lower half-plane of 
the second sheet, and the integral § • • • is clockwise around the resonance 



Sr) in the second sheet. For the ana- 



pole Sr = (Mr - tT R /2) 2 (s p 
lytic extension ( |5.2| ) to be meaningful, the integral around should van- 
ish, and the wave functions, which are determined for the physical values 
M so = [4m 2 , oo), should have a unique extension to the non-physical values 
of s, i.e., on (— oo,4m 2 ), which, as can be seen in ( |5.2| ) are always occur- 
ring in the second sheet. Since our work follows closely the derivation of the 
non-relativistic Gamow vectors, to justify (|5.2j ), we invoke below the same 
arguments used for the non-relativistic case. 

To prove the vanishing of the integral in ( |5.2| ), we make the same 
assumption about the growth of the S'-matrix as made in |25| , namely the S'- 
matrix on the second sheet is bounded by a polynomial, i.e., for large s, there 
exists a polynomial P(s) such that |S//(s)| < |P(s)|. With this assumption, 
we obtain the relation 



\ds(r(p)\j3S-)S(s)( + j 3 s\ ( j ) + (p)}\ 



< I |ds(^-(p)|j 3 S-)P(s)( + j3S|0 + (p))| 

'Coo 



(5.3) 
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From Proposition |372|, it follows that P(s)( + j 3 s|0 + (p)) G S fl H 2 _ . Hence, a 
straightforward application of Holder's inequality shows that 

^-(p)|j3S-)P(s)( + j 3 s|0 + (p)) E U\ . (5.4) 



With (|5.3| ) and (5.4), the vanishing of the integral around follows from 
Corollary |BT|. 

As mentioned above, the wave functions (pj3[sj]~|?/'~) and (pj3[sj] + |</> + ) 
should be obtained unambiguously everywhere on the real line of the second 



sheet. The Hardy class assumptions ( |3.11| ) ensure that this is actually the 



case. This follows from the remarkable property that every Hardy class 
function is completely determined from its values on a half-axis of the real 
line. In other words, there exists a bijective mapping 

e : snn 2 T ^ snn 2 T \ Rso . (5.5) 

This result, which follows from a theorem of van- Winter's theorem ( p.2| , 
Appendix B below) is as crucial here as it is for the formulation of the non- 
relativistic Gamow vectors. The 9 function in ( |5.5| ) allows the values of the 
wave functions on the non-physical region of the Mandelstam variable s on 
the second sheet to be uniquely determined from their values on the physical 
range [4m 2 , oo). 

Having elucidated how the Hardy class assumptions ( |3.11| ) provide the 



mathematical justification of (|5.2|) , we now consider the integral around the 
pole separately 



ds^-(p)\j 3 s-)S II (s)( + j 3 s\<f ) + (p)) . 



For this integral, we expand the partial ^-matrix S n (s) into a Laurent series 
about the resonance pole 

S n (s) = + A(s) , (5.6) 

s - s R 

where r is the residue of S(s) at the pole and A(s) is an analytic function. 



The pole term in (]572 ) can then be expressed as 



ds(^-(p)|j 3 S-)( + j3S|0 + (p)) 

S - S R 

= -2nir (tp- (p) \j 3 s R ) ( + j 3 s R 1 4> + (p) ) 



roc 

/ ds(^-(p)|i 3 s-)( + i3s|0 + (p)) • (5.7) 

J-ooII S — S R 
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The first equality in ( |5.7| ) the well known theorem of Cauchy and the sec- 
ond equality follows from a theorem of Titchmarsh ( (p.l| ), Appendix B). 
With (£70, becomes 

(i>is(p),4>i(p)) = -2mr (4>-(p)\j3S B )( + j 3 s R \<P + (p)) + B r M (5.8) 
where Bj 3 (p) is given by the other non-zero integrals in ( |5.2| ): 

poo 

B j3 (p) = / ds(i;-(p)\j 3 s-)( + j 3 s\ ( p + (p))S(s + te) 

J 16m 2 

fOO 

ds{r{p) I ./:;s }; './:;S | <P + (p))S n (s - ie) 

/4in J 
ds(ip-(p) | ./:;S ){+j 3 s | + (p))5 /7 (s - ie) . (5.9) 
-oo 

Since ( |5.8|) is valid for any registered two pion state if)~ + _, we can omit the 
arbitrary - G or ip~ and represent a two pion in-state density by the 
functional equation in the space $+: 



16m 2 
4m 2 



^33 

where 



[p) = -2mr\pj 3 [s R3R }-)( + j 3 s R | <f) + (p)) + \B ja (p)) , (5.10) 



poo 

\B j3 (p)) = / rfs|pj3[sj R ]-)( + j3s|0 + (p))S(s + 2e) 

J 16m 2 
/>oo 

- / rfs|pj 3 [sj«]-)( + j3S|0 + (p))^ 
J 16m 2 

/4m 2 
ds|pj 3 [sjfl]-> ( + J 3 s I <P + (P))S n (s - ie) . (5.11) 
-oo 



s — ie) 



The expansion ( [5 . 1 0| ) is the relativistic version of the complex basis expansion 
obtained for the non relativistic case (cf. equation (5.40) of 0). This is the 
complex basis expansion if there is one resonance in the partial wave. If there 
are two or iV resonances with j = j R at s = s Rl , s R2 , ■ ■ ■ (e.g., p(1450) and 
p(770) in 7r + 7r~ — > 7r + 7r~), then we obtain for each an additional integral 
around the pole in ( |5.2| ). For each of these poles separately we follow the 
above procedure and obtain, in place of the first term on the right hand side 
of ( |5.10| ) and ( |5.12[ ) below, a sum over the (superposition of interfering 
resonances) . 
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As prescribed in (|3.27| ), the in-state + of the two pions (in the subspace 
3 — 3r) is obtained from the in-state density <fi~j~ 3 (p) by integration over the 
4- velocities and summing over j'3. Thus, 

33 

= E / ^\pUsrJr]-)MP) + \B) , (5-12) 



where 
and 



;p) = -2irir( + j 3 s R \(f> + (p)) = -27rzr( + pj 3 [s R J R ]|0 + ) , (5.13) 
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The state vector Ip^s^'^] ) which appears in (|5.10 ) is the relativistic Gamow 



vector that we set out to construct. As apparent from the derivation lead- 
ing to (|5.10Q , it is obtained from the analytic extension in s of the Dirac- 
Lippmann-Schwinger kets |pj3[sj_r]~) to the resonance pole in the second 
sheet of the analytically continued 5-matrix. The first term in the right hand 
side of Q5.12p , which is a continuous linear superposition of \pjs[sRjR]~) over 
the 4- velocity p with fixed values of Sr and jn, represents a velocity wave 
packet of the resonance particle. We denote it by 



3RSR 

33 



E / Sl^3te]-)0i 3 (p)- (5-15) 



We recall that <fij 3 (p) in ( |5.15|) , defined by (|5.13|) , is a Schwartz function with 



respect to p, as required by ( p.llbj) . 



To obtain the invariant mass-square distribution of |pJ3[s_rJ,r] ), we de- 
duce from the second equality in fl5.7|) that 

w-(p)b-3s s ) = r *r*'f2y" tf)lA ° 

V ^ " " i-oo„ ( + j3S R \<t>+(P)) 5-Sh 



CO 



ds <^-(p)U3S~> (516) 
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where — oo/j signifies that the "unphysical" values of s, (—00, 4m 2 ), occur in 
the second sheet. Since (|5.16|) is valid for any registered two pion state ?P~+ n - , 
we can represent the state vector corresponding to the resonance particle by 



\P33[sr3r] ) 



i 

2^ 



ds 



-OOJJ 



s - S R 



(5.17) 



The integral representation in ( 5.17 ) shows that the relativistic Gamow vector 
has a Breit-Wigner distribution in the s-variable (invariant mass-squared) 



(pJs^Jr} iPjs^RjR] ) 



Sr 



— C&U < s < 00 



(5.18) 



Exactly as in the non-relativistic case (cf., Proposition 4 in [25]), it can be 
shown that |pj3 [shJr]~) with the integral representation ( |5.17| ) is a continu- 
ous antilinear functional on $ + , i.e., |pj 3 [s/jj/j]~) G $+. 

The relativistic Gamow vector \pj3lsRjR\~) is a generalized eigenvector 
of P M with a complex eigenvalue. To see this, we use the integral represen- 
tation ( |5.7D with the vector P^ip~ G $+: 



(PttP \P33[sr3r] ) 



1 

i 

2it 



pM$3b] ) 



—00 
00 



— / ds 



s-Sr 
V$vM~ \P3z[ s 3rY) 



S - Sr 

SrPM~ \P3z[$r3rY) ■ 



(5.19) 

In ( P-19|) we used (|3.13a|) to write {P l ^}~\Ph[sj]~) = y/sp^~\pj 3 [sj]~) 
and fl3.14p to assert that y/sp^ip^lpj^siY) is a Hardy class function from 
below, so that Titchmarsh theorem p.l| can be applied to obtain the last 
equality. Similarly, 



(M 2 ip~ I pj 3 [s R j R ] 



' 2tt 

i 

~ 2^ 



00 , (MV 
ds 



P3A s 3b\ ) 



-00 
—00 



S - Sr 
PjslsjR]') 



(5.20) 



s - s R 
P33[sr3r]~) ■ 



Equation (|5.20|) , valid for all ip G $+, is the mathematical expression that 
\P3z[ s r3r\~) is a generalized eigenvector of M 2 with the complex mass square 
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Sr as eigenvalue. This is written equivalently as 

M 2X \pj 3 [snj R ]-) = s R \pj 3 [s RjR }-) . (5.21) 

In exactly the same way, the wave-packet ( |5.15|) is a generalized eigenvector 
of the mass square operator M 2 with the eigenvalue s R : 

^4 B = ^L' (5-22) 

Hence, <fif RSR i n ( |5-15| ) represents velocity wave-packets of the unstable par- 
ticles associated with the pole Sr. 



6 Conclusion 

In this paper we have discussed a state-vector description of relativistic un- 
stable particles. Following the norms of particle physics phenomenology, we 
can anticipate at the outset that the states of an unstable particle ought 
to be characterized by the values of its mass, width and spin, the first two 
of which have been combined in several different ways to a single complex 
quantity which is associated with the position of the resonance pole of 
the relativistic S- matrix. As affirmed by Wigner's classic paper II], a sta- 
ble particle can be given a state vector description where the state vectors 
are specified by the generalized eigenvalues of momenta, real mass and spin. 
These state vectors in fact furnish an irreducible unitary representation of 
the Poincare group. 

In the theory we have proposed in this paper, the state vector description 
of unstable particles arises from the relativistic Gamow vector \pjs[sRjR]~) , 
which are vectors (associated to and) obtained from the pole term of the S- 
matrix at s = sr. It is a simultaneous generalized eigenvector of the invariant 
operators M 2 = P M P M and W = —w fJj w fl as well as the four velocity opera- 
tors P M and the third spin component W3. Evidently, the Gamow vector is a 
generalized eigenvector of the momenta as well, an immediate consequence 
of the defining identity P M = MP M = y^SRp^. Thus there exists a manifest 
parallel between Wigner's state vector description of stable particles and the 
Gamow vector educed state vector description of unstable particles. The 
fundamental difference is the complexness of the (generalized) eigenvalue of 
the invariant mass square operator M 2 in the latter case. 

8 Here w = \e^ pa P v JP a . 
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The relativistic Gamow vector |pj3[sj?jV]~) with the complex general- 
ized eigenvalue ^/sr of the (essentially) self-adjoint mass operator M was 
derived here by analytically extending the Dirac-Lippmann-Schwinger kets 
|pj3[ s Jfl]~) i n t ne Mandelstam variable s into the lower half complex plane 
from its initial range (mi + m-i) 2 < s < oo, appertaining to a scattering 
process of two particles, (m 1; ji) and (m 2 , j<i). The point Sr in the lower half 
plane at which this analytic extension yields the Gamow vector \P3z[^r3r\~) 
is where there exists a simple pole of the analytic S'-matrix designating the 
unstable particle. 

Aside from being an alluring proposition on theoretical grounds, for it 
puts our understanding of stable and unstable particles on an equal footing, 
this state vector description of unstable particles by Gamow vectors can be 
seen to unify various properties heuristically attributed to unstable particles 
and resonances. First, Eqs. ( |5.16| ) and ( |5.17| ) show that the Gamow vector 
\P3z[sr3r\~) can be resolved into an integral representation over a Breit- 
Wigner distribution — !— in the Mandelstam variable on the second sheet, 
an attribute of its association to the resonance pole at Sr. Although not 
discussed here, it will be shown in the sequel to this paper that the time evo- 
lution of the Gamow vectors is exactly exponential. Further, this exponential 
time evolution will be shown to be given by e~ VRt (with t being the time in 
the rest frame) where Sr = (Mr — ir^/2) 2 is the S'-matrix pole position. 
This implies that the width of the Breit-Wigner distribution of the particle 
is related to its lifetime t r via the exact identity Y R = we can then define 
Mr as the real mass of the resonance. Thus the characterization of unstable 
particles proposed here sheds some light on the recent debate among particle 
physics phenomenologists and theorists on the definition of the mass of the 
Z-boson and certain hadronic resonances |21 |. 



We shall derive the time evolution of Gamow vectors in the forthcoming 
sequel where we shall study the transformation properties of the Gamow vec- 
tors under Lorentz transformations and space-time translations. This study 
will show that the linear span of Gamow vectors \pjs[sRjR]~) , pet 3 , fur- 
nishes an irreducible representation of the causal Poincare semigroup, defined 
as a semi-direct product of the group of homogeneous Lorentz transforma- 
tions with the semigroup of space-time translations into the forward light 
cone. From our requirement that the complexness of the momenta be solely 
due to that of mass, i.e., by the required reality of the velocities p^ = ^==, it 
will follow that the (homogeneous) Lorentz subgroup can still be represented 
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unitarily. However, translations will only form a semigroup, not a group. 
This group theoretical study will further support our view that the charac- 
terization of stable and unstable particles need not be qualitatively different 
in that the particles of either kind can be specified according to their mass 
and spin by means of irreducible representations of the Poincare transfor- 
mations: in the stable case, we have irreducible unitary representations of 
the Poincare group, labeled my real masses; and in the unstable case, irre- 
ducible representations of the causal Poincare semigroup, labeled by complex 
masses. The semigroup of the latter case can be interpreted as representing 
the asymmetry of time evolution (say, at rest) of unstable systems and the 
irreversible character of decay processes. 
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Appendices 

A Overview of Rigged Hilbert Space Con- 
cepts 

A Rigged Hilbert Space [^] is the result of the completion of a scalar product 



space with respect to three different topologies. The completion of a vector 
space with respect to some topology r amounts to including in this space the 
limit points of all r-Cauchy sequences. If one starts with a scalar-product 
space \& and completes it with respect to the norm induced by the scalar 
product 



a Hilbert space 7i is obtained. On the other hand, if \1/ is completed with 
respect to a topology r$ defined by a countable number of norms with some 
qualifications, a countably normed space $ is obtained. This countably 
normed topology r$ is finer than the Hilbert space topology so that there 
are more r$-neighborhoods than r^-neighborhoods. Hence: 

^ C $ C U. 
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A third space of interest is the space of antilinear functionals on $, denoted 
by $ x . Since $ C H, it follows that H x C $ x . But, from Hilbert space 
theory, H = H X . Hence 

$cftc$ x . (A.l) 



The triplet QA.l ) is called a Rigged Hilbert Space when $ is nuclear and dense 



in Tt (with respect to Ty). The fundamental physical axiom of the Rigged 
Hilbert Space formulation of quantum physics is that the set of states of the 
physical system do not inhabit the entire Hilbert space Ti but an appropri- 
ately defined dense subspace <3> of Ti. The countably normed topology of 
<3> is constructed so as to yield the algebra of relevant physical observables 
continuous as mappings on $. It is this feature of Rigged Hilbert Space 



theory that is made use of in Section |3J| in making the distinction between 



the set of prepared states and registered observables $+ by taking §± as 
dense subspaces of the same Hilbert space Ti as in (|3.6| ). This distinction is 
what allows semigroup time evolution to be incorporated into the quantum 
mechanical theory. 

The action of an element F G $ x on G $, F(<p), is denoted-in the 
Dirac bra-ket notation-by 

F{<f>) = (<f>\F) . 

Since Ti C $ x , it follows that the Dirac bra-ket (<p\F) is an extension of the 
Hilbert space scalar product in the sense that 

= for FeH. 

The topology on <3> x , denoted by r$x, is the weak* -topology induced by $ 
on $ x . This means that convergence in $ x is defined by 



F ^ ((t>\F) -> (4>\F) , for all <p G $ . (A.2) 

To every r$-continuous operator A on there corresponds a r^x -continuous 
operator A x defined on $ x by 

(0|A x F) = (A(j)\F) , for all G $ , F G $ x . (A.3) 

The operator A x is called the conjugate operator of A. It is an extension of 
the Hilbert space adjoint operator A\ since for F G Ti we have 

{(f)\A x F) = (A(j),F) = (4>,A^F) for FeH. (A.4) 
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Hence, 



A% cA^cA x . (A.5) 

It should be stressed that the conjugate operator A x can be defined as a 
r$x -continuous operator only when A is a continuous linear operator on $. 
In quantum mechanics, it is impossible (empirically) to restrict oneself to 
continuous (and therefore bounded) operators A in TC. However, one can 
restrict oneself to algebras of observables {A, B, ■ ■ - } described by continuous 
operators in $, if the topology of $ is suitably chosen. Then, A x , B x , ■ • • 
are defined and continuous in $ x . 

A generalized eigenvector \F) of a 7$- continuous operator A with a gen- 
eralized eigenvalue cu G C is defined by the relation 

(A(j)\F) = (<j)\A x F) = u(<t>\F) , for all G $ . (A.6) 

Since the vector in ( |A.6|) is arbitrary, ( |A.6| ) can be formally expressed as 



A x \F)=to\F). (A.7) 
In the Dirac notation the x in ( |A.7p is suppressed so that ( |A.7p reads 



A\F)=u\F). (A.8) 

If A is a self-adjoint operator, suppressing the x as in ( |A.8| ) does not lead to 
confusion since A = A* C $ x . However, if A is not self-adjoint, a clear dis- 
tinction between the operator and its conjugate should be made. The concept 
of generalized eigenvectors ( |A.7| ) in Rigged Hilbert Space mathematics allows 



the description of "eigenstates" which do not exist in the Hilbert space. For 
instance, the Dirac scattering kets are generalized eigenvectors with eigen- 
values belonging to the continuous spectrum, and they are not Hilbert space 
elements. The Gamow vectors, which are used to describe decaying states, 
are also generalized eigenvectors which are not in 7i, but, unlike in the case 
of scattering states, their complex eigenvalues do not belong to the Hilbert 
space spectrum of the Hamiltonian. 



B Hardy Class Functions on a Half-plane 
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Definition B.l (H.± 1 < p < oo). A complex function f(x+iy) analytic in 
the open lower half complex plane (y < 0) is said to be a Hardy |27j class 
function from below of order p, 7i p _, if f(x+iy) is L p -integrable as a function 
of x for any y < and 



sup / dx \f(x + iy)\ p < oo . (B.la) 
y<0 J-oo 

Similarly, a complex function f(x + iy) analytic in the open upper half com- 
plex plane (y > 0) is said to be a Hardy class function from above of order 
p, Ti p +1 if f(x + iy) is L p -integrable as a function of x for all y > 0, and 



sup / dx \f(x + iy)\ p < oo . (B.lb) 
y>o J-oo 

A property of 7i p ± functions is that their boundary values on the real axis 
exist almost everywhere and define an L p -integrable function, i.e., if / G 7~t±, 
then its boundary values f(x) G L P {R). Conversely, the values of any H± 
function on the upper /lower half-plane are determined from its boundary 
values on the real axis. This result is provided by a theorem of Titchmarsh: 

Theorem B.l (Titchmarsh theorem). If / G H^, then 

-1 f°° f(t) 

/(*) = T^/ r 1 ^' forlmz<0, 
2m ./_oo t - z 



and 

/(*) 



oo 



t - Z 



dt = , for Im z > 



Similarly, if / G H+, then 



i r°° f(t) 

f(z) = / ^±±dt, forlm^>0 

2m ./_oo * - z 



and 

/(*) 



dt = , for Im z < . 



This one-to-one correspondence between the ?{± functions and their bound- 
ary values on K. allows the identification of f(z) with f(x) for / G 1-L P ±. 

The following results are related to the decay properties of the Hardy class 
functions. They are straightforward generalizations of the corresponding 
results of [^| and are needed for the construction of the relativistic Gamow 
vectors. 
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Proposition B.l. Let be the infinite semi-circle in the lower half com- 
plex plane. If / G UlL, then 



L 




= 




z 





Proof. Let C r be the arc with radius r shown in Figure |3|. Then 



-dz 



< I \f (re w )\de = ^ > r \f(-re w )\d6. 

C r Jl/r 



Since / G Ti^, then there exists C such that 



|/(-re")| < 



C 



Thus 



/(*) 



dz 



< 



r sin 0)Vp 
2C r /2 1 



(cf. p|] page 149) 



-l/p 



l/r 



(sin^Vp 



d0. 



Using 



sin# > - 9 3 /Q > 9(1 - tt 2 /24) , for l/r < 9 < tt/2 
we obtain for (|B.2j) 







< 




z 





2C 



■"' 2 d9 



r i/ P (i _ vr 2/24) 1 /p y 1/r 0i/p 

_ 2C J log (f ) 

~ (1 - 7 r 2 /24) 1 /p r 1 /p | ^ [(f) 1 " 1 /* - (i) 1_1/p ] 

Therefore, as r — > oo, we obtain 

0. 



(B.2) 



p — 1 

1 < p < oo 



(B-3) 



/joe 






z 



□ 



Corollary B.l. Let / G 5n W 2 , 5 G 5nW 2 , then 

|/(z)<7(z)dz|=0. 
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Proof. Since / G S DT-Ll, then xf(x) e 5flW! @. A straightforward 
application of Holder's inequality shows that xf(x)g(x) G H}_. Then, from 
the above lemma 



\f(z)g(z)dz\ 



zf{z)g{z) 



dz 



0. 



□ 



A remarkable property of Hardy class functions that is used in |25) is that 
they are uniquely determined from their boundary values on a semi-axis on 
the real line. This result is provided by a theorem of van Winter |29| . Before 



stating the van Winter's theorem below, we define first the Mellin transform 



Definition B.2 (Mellin transform). Let f(x) be a function on 
Mellin transform is a function defined almost everywhere on M as 



.. Its 



His) 



f(x)x is - 1/2 dx. 



W 2 Jo 

provided that the integral exists for almost all s G 



Theorem B.2 (van Winter). A function f(x) G L 2 (R + ) can be extended 
to M~ = (—00, 0] to become a function in Ti 2 ^ if and only if its Mellin trans- 
form satisfies 

/oo 
(l + e 27TS )\H(s)\ 2 ds< 00. 
-00 

This extension is unique. The values of f(z) for z = pe ld for < 6 < n, 
p > are given by 

In particular for negative values of x, f(x) is given by 

00 

H(s)(xe t7T )- ts - 1/2 ds. 

00 

A similar result can be obtained for 7i 2 __ . 
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Figure 1: Analytic properties of S(s) 
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Figure 2: Contour for the analytic continuation of ( |3.28| ). 
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